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1 Introduction 

Determinantal (a.k.a. fermion) random point processes were introduced in 
probability theory by Macchi about thirty years ago (QSIj El)- In the 
last ten years the subject has attracted a considerable attention due to its 
rich connections to Random Matrix Theory, Combinatorics, Representation 
Theory, Random Growth Models, Number Theory and several other areas of 
mathematics. We refer the reader to the recent surveys ([201; 0) |H1)) and 
research papers on the subject {[T], 0, 0, 0, [El, t^, Hi, [III, P| . 

In this paper we shall consider determinantal random point processes on 
the real line with the translation-invariant correlation kernel. In other words, 
a one particle space X is given a.s X — Ti}, and the space of elementary 
outcomes f2 consists of the countable, locally finite particle configurations on 
the real line 

n = = {x,)+^_^ : #{x,e[~N,N])<+^, ViV>0}, 

where x, G 7^^ i = 0,±1,±2,..., and #(a;, £ [~N,N]) =: #([-iV,iV]) 
denotes the number of the particles in the interval [— A^, A^]. Let us denote 
the set of the non-negative integers by = {0, 1,2, . . .}, and the set of 
the positive integers byA/'={l,2,...}. The cr-algebra J- of the measurable 
subsets of n is generated by the cylinder sets C^^f^'"'^^ = : — 
rij, = 1, . . . , fc}, where k is an arbitrary positive integer, k G A/', /i, . . . , 
are arbitrary disjoint subintervals of the real line, and ni,n2, ■ ■ ■ ,nk G Z^. 

A probability measure V on the measurable space {f2, J-) defines a random 
point process {Q,J-,V). A random point process is called determinantal if 
its fc-point correlation functions have determinantal form 
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Pk{xi, . . . ,Xk) = det{K{xi,Xj))ij=i^,,,^k, A: = 1,2,..., (1) 

where K{x,y) is usually called the correlation kernel of the determinantal 
random point process. We remind the reader that /c-point correlation functions 
are defined in such a way that 

eY\_#{Ii) = pk{xi, . . . ,Xk)dxi . . .dxk (2) 

for the arbitrary disjoint intervals /i, . . . ,Ik- Since the r.h.s. of Q is non- 
negative, it follows that the correlation kernel K{x,y) has non-negative mi- 
nors. If, in addition, the integral operator K : L'^iji}) L'^{TZ^), {K f){x) — 
K{x,y)f{y)dy, is Hermitian, one can conclude that K is non-negative 
definite, i.e. Spec{K) G [0, +oo). In the Hermitian case one can show that 
the necessary and sufficient condition on K to define a determinantal random 
point field ^ is 

< < 1, (3) 
in other words both K and 1 — K must be non-negative definite operators 

(im, csi). 

In this paper we consider the translation-invariant kernel 

r+oo 

K{x,y) = g{y — x), where g{x) ^ I eyip{2Trixt)(l){t)dt, (4) 



and (l)(t) is an even real-valued integrable function 

cj,eL\R'). (5) 

It follows from Q that 

< 4){t) < 1 (a.e.). (6) 
In addition, we assume that the following technical conditions are satisfied 

+2 



t^(l){t)dt < +00, (7) 
\9{x)\ < (8) 

1 -I- |X| 2+*^ 

where C is a positive constant, and e is an arbitrary small positive constant. 

Let 1/ be a large positive number. Consider a restriction of a configuration ^ 
to the interval [0, L] . Let us denote the points of [0, L] by Xi , X2 , . . . , ai^rji) , 
where N{L) is the cardinality of f n [0, L]. We assume that the points in 
^ n [0,L] are ordered: xi < X2 < ■ ■ ■ < a:^Ar(L)- It is a well known (see 
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e.g. 1201) that with probabihty 1 no two particles of a determinantal random 
point process coincide. We are interested to study the nearest spacings 6i — 
Xi^i — Xi, i = 1, . . . , N(L) — 1, between the neighboring particles. Functional 
Central Limit Theorem for the empirical distribution function of the nearest 
spacings of particles in [0, L] (in the limit L oo) was proven in |21| for 
K{x,y) — (i.e. is the indicator of [— i,i]), and for similar kernels 

arising in Random Matrix Theory. It was shown in 20 that the result could be 
extended to a quite general class of translation-invariant correlation kernels. 

In this paper we study the smallest nearest spacings. Our main result is 
the following 

Theorem 1. Let {f2,J-,'P) he a determinantal random point process on the 
real line with the translation- invariant correlation kernel K{x, y) = g{y — x) 
satisfying ^-0). Then the number of the nearest spacings less than s/L^^^ 
in the interval [0, L] converges in distribution to the Poisson random variable 
with the mean as'^, in the limit L — > oo, where 

a = xg(0).9"(0) = — / mdt / t^mdt. (10) 

Let r]{L) = L^/^ min{9i, i — 1, . . . , N{L) — 1}. In other words, ?7(i) is the 
smallest nearest spacing in [0,i], rescaled by L^/"^. Theorem 1 immediately 
implies 

Theorem 2. Let the conditions in Theorem 1 he satisfied. Then 

lim Pr(7/(L) > s) = exp(-as^). (11) 

L — ^oo 

The method of the proof of Theorems 1 and 2 relies on the detailed analysis 
of fc-point correlation and cluster functions of the s-modified random point 
process, introduced in [201, [21] • 

The rest of the paper is organized as follows. Point correlation and cluster 
functions, and s-modified random point processes are discussed in Section 2. 
The proofs of Theorems 1 and 2 are given in Section 3. 

We will use the notations const, constk, Const, to denote various positive 
constants throughout this text. The values of these constants may be different 
in various parts of the paper. We shall use the notation / = 0{g) if the ratio 
f/g is bounded from above and below by some positive constants, and the 
notation / = o{g) if the ratio f/g goes to zero. 

2 Correlation and Cluster Functions 

We start this section by recalling the definition of a fc-point cluster function 
(sometimes also known as the Ursell factor). For additional information we 
refer the reader to [T?)j. [TT]. Pj. 
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Definition. The l-point cluster function ri{xi, . . . , xi), Z = 1, 2, . . . , of a ran- 
dom point field is defined in terms of the point correlation functions by the 
formula 

m 

n(xi, . . . , Xi) = ^(-l)™-i(m - 1)! n P\gMG^)) (12) 

G j=l 

where the sum is over all partitions G of [I] = {1,2,...,^} into subsets 
Gi, . . . , Gm, m = 1,. . . ,L and x{Gj) = {xi : i e Gj}, \Gj\ = #(Gj). 

The point correlation functions can be expressed in terms of the point 
cluster functions as 

■m 

p,(xi,...,xO = ^nnG,i(s(G,))- (13) 

G j=l 

The reader can observe that (|12|l is the Mobius inversion formula applied to 
(|13|) . In particular, 

pi{x) = ri{x), 

P2ixi,X2) = r2{xi,X2) + ri{xi)ri{x2), 

P3{xi,X2,X3) = r3{xi,X2,X3) + r2 (xi , X2 )ri (xa ) + r2(a:i, a;3)ri(x2) + 
r2{x2,X3)ri{xi) + ri{xi)ri{x2)ri{x3). 

It follows from (|13|) and ^ that for determinantal random point fields 

ri{xi, . . . ,Xl) = {-ly^^ ^ K{xi,X„(^l))K{x2,Xa(2))---K{xi,Xa(l)), 

cyclic aeSi 

(14) 

where the sum in (|14|l is over all cyclic permutations. In other words, for 
determinantal random processes the difference between the formula (|14|l for 
the Z-point cluster function and the formula 

Pl{xi, ...,Xl)^ ^ {-iyK{xi,Xa(l)) ■ K{x2,Xa[2)) ' • • ■ ' K{xuXa(l)) (15) 

for the Z-point correlation function is that in (|15|l the summation is taken 
over all permutations in Si, and in 114() the summation is over the cyclic 
permutations only. One can rewrite 114(1 as 

n{xi,...,xi) = (-l)'"^-y X! -^(^<7(i)'2;£7(2))^(2;fT(2),a;„(3))...^(a;<T(i),2;CT(i)). 

(16) 

It follows from ^ that the integral of the /c-point correlation function 
over the /c-dimensional cube [0, L]'^ is equal to the fc-th factorial moment of 
the counting random variable #(/), / — [0, L], namely 
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E#{I)mi) - 1) . . . (#(/) -k + l)= f pk{x,, . . . , Xk)dx, . . . dxk. (17) 

The integral of the fc-point cluster function, in turn, can be expressed as a 
linear combination of the cumulants of #([0, L]). Namely, let Cfc(L) denote 
the fc-th cumulant of #([0, L]), and Vfc(L) — J^^ rk{xi,X2, ■ ■ ■ Xk)dxi x 
dx2 ■ ■ ■ dxk- Then (see e.g. 0, |23) 

n— 1 n— 1 

To apply the machinery of the cluster functions to the problem at hand, 
we consider a so-called s-modified random point process, which can be con- 
structed in the following way. We start with a random configuration ^ = 
{^i)t^oo fronr the original random point field, and keep only those points 
Xi for which there is exactly one neighbor to the right within distance s, i.e. 
Xi+i — Xi < s, Xi+2 — Xi > s. The points Xi for which this conditions is 
not satisfied are thrown away. As a result, we obtain a new random config- 
uration ^(s) C ^, such that ^(s) = {xi : Xi+i — Xi < s, Xi+2 — Xi > s}, 
where . . . < x^2 < X-i < xq < xi < X2 < ■ ■ ■ are the points of the original 
configuration ^. The number of spacings less than s of the original random 
point field in the interval [0, L] is related to the number of points of the s- 
modified random point field in [0, L]. As will be shown later, for large L and 
s ^ L^^/'^, these two counting random variables coincide with probability 
very close to 1. 

Since the moments and the cumulants of the counting random variable 
#([0, L]) can be expressed in terms of the integrals of point correlation and 
cluster functions H17(l . H18(l . it is essential to be able to calculate the point 
correlation and cluster functions of the s-modified random point process. 
We shall denote the k- point correlation and fc-point cluster functions of 
the modified random process by pk{xi, . ..,Xk',s) and rk{xi, ... ,Xk',s), cor- 
respondingly. It follows from the inclusion-exclusion principle that, provided 
\xi — Xj\ > s, ^ < i ^ j < k, one obtains 



, V- (-1)" r ' f f 

Pk{xi,...,xk;s)= ... / 

„j_Q Jxi Jxk J I(xi,...,Xk\s)™- 

P2fe+m(a;i, . . . , Xfc, 2/1, . . . , j/fc, zi, . . . , z„r)dyi . . . dykdzi . . . dzm, (19) 

where /(xi, . ..,Xk;s) = U™ J^^i, Xi+s], and I{xi, ... ,Xk; s)™ = /(xi, . . . s)x 
... X I(xi , . . . , Xfc ; s) stands for the m-th fold Cartesian product of /(xi , . . . , Xfe ; s) 
(see e.g. ED). 

In the determinantal case the formula for the fc-point cluster function 
of the s-modified random process has a somewhat similar structure ( 20 , 21 ). 
Provided \xi — Xj\ > s, 1 < i ^ j < k, one has 
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rkixi, ...,Xk;s)= V 



m— 



JKfc J I{xi,...,Xk;s) 



P2k+mixi, • ■ • , a;fc, yi, • ■ • , 2/fc, zi, . . . , Zm)d2/i . . . dykdzi . . . dzm, (20) 

where 

P2k+mixi, ■ ■ ■ ,Xk,yi, ■ ■ ■ ,yk, Zl, . . . , Zm) = ^ {-iyK{xi,a{xi)) . . .K{xk,<7{xk)) X 

''■£5'2A: + m 

K{yi,a{yi)) . . . K{yk, a{yk))K{zi, a{zi)) . . . K(z^, a(z„0), (21) 

where the summation in H21|) is over the permutations a G S2k+m satisfying 
the property A described below (we note that cr acts on the set of 2fc + rn 
variables xi, . . .,Xk,yi, . . ■,yk,zi, ■ ■ --Zm)- 
Property A 

Consider an index I < i < k. Define X{i) to be the subset of the set of 
variables {xi, ... ,Xk,yi, ■■■ ,yk, Zi, ■■■ Zrn}, that consists of Xi, yi, and those 
of the variables Zi, . . . , Zm, that belong to [xi, Xi+s] (we remind the reader that 
each zi, 1 < I < m, belongs to exactly one interval [xj,Xj + s], 1 < i < fc)- 
Then for any pair of disjoint integers 1 < i ^ j < k there exists a positive 
integer r = r{i,j) > 0, such that a''{X{i)) f] X{j) ^ 0. 

We would like to bring to the reader's attention the fact that the relation 
between (I19|l and H2U|I is, in a sense, quite similar to the relation between H15|) 
and (dH). 



3 Proof of the Main Result 

The strategy of the proof is the following. We consider the rescaling s = s x 
L~3, (we shall show that the smallest spacings in the interval [0, L] are of 
order L~^/^). We shall keep s fixed as L ^ oo, so s will be proportional to 
L~3 . We are interested in the asymptotics of the integrals 



Vfe(L) = / rk{xi,X2, . . .Xk;s)dxidx2 . . .dxk. (22) 

We claim that limL^ooVi(L) = as'^, where a is defined in and 
limi^oo Vk{L) = 0, for fc > 1. 

Lemma 1. Let Vk{L) be defined as in \2'l\j , where rk{xi, X2, ■ . . Xk', s) is the k- 
point cluster function of the s-modified random point process introduced above 
and s = sL^3 . Then 

T Tr/T\ ( as^ if k = 1, 



where a has been defined in flUl) . 
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The result of Lemma 1, combined with H18|) . imphes that the number of 
the points of the s- modified random process in the interval [0, L] converges 
in distribution to the Poisson law as L ^ cxd. 

Once Lemma 1 is proven , we shall show that the number of points in [0, L] 
of the original determinantal process that have at least two neighbors to the 
right within distance s/L}^^ is zero with probability very close to 1, provided 
that L is large and s stays finite. 

Proof of Lemma 

We start with Vi{L). Consider the one-point correlation function (inten- 
sity) of the s-modified point process pi{x\ s). 



(2;;s)=V' / •■•/ prn+2ix,y,zi,...,z„,)dyd"'z. (24) 



Pi 

We claim that in the determinantal case 

px+s 

pi{x;s)= p2ix,y)D{x,y;s)dy, (25) 

J x 

where D{x, y; s) is the Fredholm determinant of the integral operator K on 
L'^{[x,x + s]), 

D{x,y; s) = det{l - K), K : L^{[x,x + S]) ^ L^{[x,x + S]). (26) 

The kernel K{u,v) in depends on x and y, and is given by the 

formula 

K{u^ v) — K{u, v) — K{u, x)Tii{x, y)K{x, v) — K{u, x)Ti2{x, y)K{y, v) 

~ Kiu, y)T2i{x, y)K{x, v) - K{u, y)T22{x, y)K{y, v), (27) 

where 

'Tii{x,y) Ti2{x,y) \ _ ( K{x,x) K{x,y)^ 



T2i{x,y)T22{x,y) J \K{y,x) K{y,y) 

Indeed, let us introduce the notation K[xi, . . . , Xk] ■— dct{K{xi, Xj))ij=i^,,,^k- 
Then, 

K[x, y, zi, . . . , Zm] = K[x, y]K[zi, z,„] = p2{x, y)K[zi, . . . , z„]. 

In other words, the conditional distribution of a determinantal random point 
process with the correlation kernel K, given there are two particles at x and 
y is again a determinantal random point process (on TZ^ \ {x, y}) with the 
kernel K (see e.g. JHl)- This allows us to rewrite H24|) as 



px+s +00 f _^\m p 

pi{x-~s)= p2ix,y){J2-^ K{zi,...,zm)d"'z)dy, (28) 

Jx J[x,x+s]"^ 
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and follows. 

The intensity pi{x;s) is constant (i.e. it does not depend on x) in the 
translation-invariant case. To estimate pi{x; s) — pi{0; s), we note that 

Pm+2ix, y, zi, . . . , z„i) = K[x, y,zi,. . Zm] < K{x,x)K{y, y)K{zi,zi) . . . K{zm, Zm) 
<g{or+\ (29) 

since the determinant of a non-negative definite matrix is bounded from above 
by the product of the diagonal entries (for the generalization of this result see 
Lemma 2 below). Then 



/•X + S fX + S px + s 

Pi{x;s)^ P2{x,y)dy~ / pj,(x,y, zi)dydzi 

J X J X J X 

Y rx-\-s nx-i-s px-^s 

-J J J P4{x,y,zi,Z2)dydzidz2 + 

^ / •■•/ pm+2{x,y,zi, . . . ,z„i)dyd™'z = 

m=3 "'^ "'^ 

CX-\-s /'X-\-s px-'rs 

P2{x,y)dy- / p3{x,y,zi)dydzi + 

J x J X 
fX + S i*x + s i*x + s 



Y px-\-s px-ts px-^s 
2 Jx Jx Jx 



{x,y, zi, Z2)dydzidz2 + 0{s ). 



To estimate p3{x,y,zi)=K[x,y,zi] and p4ix,y, zi, Z2) = K[x,y, zi, Z2] 
we recall that the point correlation functions are given by the determinants, 
and subtract the first column from the other columns both in K[x,y,zi] and 
K[x, y, Zi, Z2\- Since y & {x,x + s\, Zi & [x,x + s], i>\, and the first deriva- 
tive of g is uniformly bounded, we observe that p3(x,y,zi) — K[x,y,zi\ — 
O(s^), p4(x, y, zi, Z2) = -fC[x, y, zi, Z2] O(s^), and, therefore 

Pi{x-~s) = P2ix,y)dy+0{S^) = f\g\0)-g^{t))dt+O{S^) ^ as^'+OiS^), 



(30) 

where a has been defined in H1U|) . It follows from (|30|l that lim^^cx; ViiL) = 
Uuil^oqJq Pi{x\s)dx = lim^^oo Pi(0; s)L = as^ . 

Next, we show that lim^^oo T4(L) =0 for fc > 1. We remind the 
reader that Vk{L) has been defined as Vk{L) = rk{xi, X2-, ■ ■ ■ Xk] s) x 

dxidx2 ■ ■ ■ dxk- We start with the case k = 2. Recall (see (|2U|l ) that for 



ixi,X2,yi,y2,Zi,...,Zm) 



J-fZ f-^T]'"^ fXl+S PX2+S /■ 

r2{xuX2;S)^Y.-^ / / PT+n 

m=0 •'^i •'^2 J I{xi,X2;s)™- 

dyidy2dzi . . . dzra, (31) 
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where P4™m(a;i, a;2, yi, y2j ^i: ■ ■ • j ^m) has been defined in H21|) . 

As described in the Property A (right after the formula (|21|l '). in or- 
der to define P4™^(a;i, 2:2, yi, j/2, -^i, . . . , z^) one introduces a partition 
X{1)\JX{2) ^ {xi,X2,yi,y2,zi,...,Zm}, where X{1) consists of xi,yi, 
and those of the variables zi,...,Zm that belong to [a;i,a;i + s], and 
X{2) consists of X2,y2, and those of the variables zi,...,Zm that 
belong to [0:2, a;2 + s]. Let X{l)f]{zi,...,Zm.} = {zi^, ■ ■ ■ , z^,}, and 
X{2)f]{zi, . . .,Zra} = {z-i,,. . . Then 

pT+^{xi,X2,yi,y2, zi,..., 2,„) = K[xi,yi, X2,y2,zi, . . . , z^] - 
K[xi,yi,z^^,..., z,^]K[x2, y2,Zj,,..., Zj^_,]. (32) 

We claim that 

|pr+™(^i,^2,yi,2/2,^i,...,^™)| <(m + 4)!(CS)2+™— — (33) 

1 + |Xl — X2\^ 

where const is a constant that may depend on s, and C is the constant intro- 
duced after the formulas ©, (HJ- 

The factor (Cs)^^™ in H32() follows from the uniform bound on the 
derivative of and the fact that the m+2 variables yi, 2/21 i • • ■ 7 are within 
distance s from either xi or X2- In other words, one can subtract the first 
column in the matrices in K[xi, yi, Zi-^ , . . . , Zi^] and K[x2, j/2, ^ji , ■ • • , 
from the other columns, and subtract the first and the third column in the 
matrix in K[xi,yi, X2,y2, zi, . . . , Zm] from the corresponding columns. Such 
linear operations do not change the value of the determinants, and the new 
matrices will contain the terms g{u — w) — g{xj — w), in all columns, except 
those corresponding to xi and X2, where j = 1, 2, and u G [xj, Xj + s]. Such 
terms can be estimated from above by (max^g[^^_j„ \g'{x)\)s. 

It follows from the definition that p*'l^^{xi,X2, yi,y2, zi, . . . , Zm) can be 
written as a sum over at most (m -I- 4)! permutations, each term being a 
product m + 4 factors. As we just showed, m + 2 out of those m + 4 factors 
can be estimated in absolute value by Cs. Moreover, Property A implies that 
at least two factors in each term must be given either by g{xi — v), or by 
g{u — v) — g{xi — v), or by g{x2 — u), or by g{v — u) — g{x2 — u), where 
u £ [xi,xi + s], V G [x2,X2 + s]. The inequalities (jSJ , ijHl imply that these 
two factors each contribute an upper bound — 1 — and the desired 

l+(|a:i-a;2|-S)3+' 

estimate follows. 

We recall that we defined above /(a:i,a;2;s) = [a:i,a;i + s] (J [^2, 2:2 + s]. 
Then 



, Zm)dyidy2dzi . . .dzm\ < 
(34) 



fT"iJT} { 

P4+,n{Xl,X2,yi,y2,Zi, 

I I(xi,X2\s)" 

(m + 4)!(Cg)^+"'(2g)^+"' ^""^^ 

l + \xi- X2V+^ 
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and 



|r2(xi,X2;S)| < g -l(m + 4)!(Cgf+"'(2gf+"' ,^|""'' < 
^-^ mi 1 + \xi — ^2^^ 

m—O 



1 + \xi - X2 



1+e ■ 



We remind the reader that l|35(l has been derived for \xi — X2\ > s. 
Since s = sL"^^^, it foUows from the above estimate that 

L rL 

^4 ) 



hm / / r2{xi,X2;s)xD{xi,X2)dxidx2 ^ lim 0(s L) = 0, (36) 

where D = {{xi,X2) ■ |a;i — a;2| > sL~^^^}. 

To estimate ^2(2:1, X2', s) on D'^, we note that 

r2{xi,X2;s) = P2{xi,x2;s) - pi{xi; s)pi{x2; s). 

It follows then from (|Sn|l that pi{xi;s) = pi{x2;s) < consts^L~^. To 
estimate ^2(2^1, 2:2; s) we can assume without loss of generality that X2 < 
Xi < X2 + s. Then 

P2(a;i,X2; s) = ^ —j / /93+„(a:i, a;2, yi, zi, . . . , z^) 

dyidzi . . .dzm, (37) 

where /(a::i,X2;s) = [a;i,a;i + s]lj[^2,a;2 + s] = [x2,xi + s]. Subtracting the 
first column in p3+m{xi,X2,yi, zi, . . . , Zm) = K[xi, X2,yi, zi, . . . , Zm] from 
the other columns and using 0, we see that p3+m{xi, X2,yi, zi, . . . , Zm) < 
(m + 3)l{const * 5)™+^. Integrating over yi, zi, . . . , and summing over m 
we obtain 

P2ixi, X2', s) < consts^, (38) 

which implies 



L nL 



lim / / r2{xi,X2;s)xD'={xi,X2)dxidx2 



lim / / p2{xi,X2;s)xD'={xi,X2)dxidx2 - 

lim / / pi(xi; s)pi(2:i; s)x£,c (xi, a;2)dxidx2 = 
L^°°Jo Jo 

lim (0(S*)L - 0{f )L) = 0. (39) 

L—>oo 

Combining (j36|l and (j39|) one obtaines lini/^^oo ^2(^) = 0- 
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The argument in the case of general fc > 2 is quite similar. Again, we 
first estimate rf^{xi,X2, ■ ■ ■ ,Xf:;s) on D ^ {{xi,X2, ■ . ■ ^x^) : \xi — Xj\ > 
s, 1 < « 7^ J < fc}- We will use formulas (|20|) and H21|) . To estimate 
P2fcTm(^i' ■■■,Xk,yi,.--,yk,zi,..., Zm) we consider the partition X{1) |J • ■ ■ U ^(fc) 
{xi, . . . ,Xfc,yi, . . ■,yk,zi, . . .,Zm}, where X{i) = 

{xi, . . . ,Xk,yi, ■ ■ ■ ,yk, zi, . . . , Zm} f][xi,Xi + s]. Let X{i) f]{zi, . . . , Zm} = 
. . . Then X{{) = {x,,y,,z['\ . . . 

It follows from Property A and the inclusion-exclusion principle that 

piT+rni^i, ...,xk,yi,...,yk,zi,..., Zm) = Y.^-iy{] - 1)! n Ki, (40) 

G 1=1 

where the summation is over all partitions G = Gi |J . . . |J Gj of [fc] = 
{l,2,...,fc}, j = l,...,fc, and Ki ^ K[xi,yi, z[^\ . . . , zi^^ : i e G;]; in 
other words, Ki depends on the variables from Uj^Q^ A'(i), and it is given 
by the determinant of the matrix built from the correlation kernel K{x,y). 
We claim that 

\P2k+m{xi,...,Xk,yi,...,yk.zi,...,Zm)\ < (m + 2fc) ! (GS)'^+™ X 

/ const const const ^ , , , 

( 1 1 1 \- ■■■ ), (41) 

'1 + -a:2|5+^ 1 + 1x2 -2:3|2+^ I + \xk - xi\2+'' 

where the summation in the last factor of the r.h.s. of H41() is over all (fc — 
1)! cyclic permutations (for example, the first term in the sum corresponds 
to the cyclic permutation 1-^2^3^. ..^fc^l). We claim that 
the estimate (|1T|| follows from (jHJ, ©, ipn|l and Property A. As in the case 
fc = 2 discussed above, we use the fact that each of the fc + m variables 
2/1, ... , yk, zi, . . . , Zm lies within distance s from one of the cc^'s, i = 1, . . . , fc. 
In each Ki = K[xi, yi,Zi \ . . . , Zn} : i G G;] in l|in|l we subtract for each i G Gi 
the column corresponding to Xi from the column corresponding to yi and from 
the other columns corresponding to the variables from X{i). These linear 
operations do not change the values of determinants, and, therefore, do not 
change the value of P2fcTm(^i' ■ ■ ■ ^ ^k,yi, ■ ■ ■ ,yk, zi, . . . , z,„). Now, according 
to the Property A, we observe that P2fc+m ^ '^^ most (TO + 2fc)! terms. 
Each term is a product of m + 2fc factors. Property A assures that each term 
in the sum can be put into correspondence with a cyclic permutation a on the 
set of fc variables xi, X2, . . . , Xk, in such a way that fc out of m-|-2fc terms in the 
product are of the form g(a;(j(i) — u), or g{u — v) — g{xcr[i) — v), i = l,...,fc, 
where v G [a;<j(i+i), a;cr(i+i) + s], and a{k -f 1) = fT(l). The bounds (O, © 
then imply (|41|l in the same manner as has been shown in the case fc = 2. 
Therefore, 

■■■ / plT+mi^i: ■ • • , yi, . . . , zi, . . . , Zm)dyi ...dyk 

dzi... dzm\ < (to + 2fc)!(G5)'=+™(fcs)'=+™ x 
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, const const const 



^l + \xi-X2\i+'l + \x2-X3\i+' 1 + \xk - xi\i+'' 

provided x = (xi, . . . ,2;^,) G D, i.e. \xi — Xj\ > s for i ^ j. Then on D we 
have an estimate 



(42) 



+00 

|rfc(xi,...,xfc;5)| < V — (2fc + m)!(Cs)'=+™(fcS)'=+™ 



const const const 



X 



m! 

m— 



< 



1 + \xi — X2\^'^'^ 1 + 1^2 — a^sl 2+^ 1 + jxfc — a:i| 2+^ 

^ / ^9l./ const const const , 
Constk{s) { 1 — X J — X J h ...), 

^l+\Xl - X2\^+' 1+|X2-X3|2+' 1+ -_a;^|2+^ 

and 

rk{xi, . . .,Xk;s)xDixi, . . . ,Xfc)dxi . ..dxk\ < Constk{L^+'^''-'^'>^i-'^)s^'' L^^ 

(43) 







The last estimate implies 



L nL 



lim / .../ rk{xi,...,Xk;s)xD{xi,...,Xk)dxi...dxk = 0, (44) 



L^o^jQ JO 

for fc > 3. Our next goal is to show that 

L 



lim / .../ rk{xi, . . . ,Xk;s)x'bixi^ ■ ■ ■ :Xk)dxi . . .dxk = 0. (45) 

To estimate rk{xi, . . . , Xk', s) on D^, we rewrite the formula (|12|l that expresses 
the fc-point cluster function in terms of point correlation functions: 

rk{xi, ...,Xk;s)= pk{xi, . . . ,Xfc; s) - pi{xi;s)pk-iix2,X3, ...,Xk;s)- 
Pi{x2] s)pk-i{xi,X3, . . . ,Xk;s) . . . - pi{xk;s)pk-i{xi,X2, ■ ■ .,Xk-i;s) + 

2p2{xi,X2; s)pk^2{x3, ...,Xk;s)+ 2p2{xi,X3; s)pk^2{x2,X4, ...,Xk;s) + 

. . . p2{xk-i,Xk\ s)pk-2{xi,X2, ■ . . ,a;fe_2; s) ~ ... (46) 

We claim that the integral of each of the terms in (|46|l over [0, L]™ p| D'' has 
a zero limit as L — > 00. To prove it, we consider an arbitrary term in (|46|1 . 

PkAxi,X2, . . . ,Xk^,s)pk.Axk^+i,Xk^+2, ■ ■ .,Xk^+k2\s) ■■■ Pki (a;fei+...fe,_i+i, . . .,Xk]s), 

(47) 

where ki + k2 . . . + ki — k, ki > 1, i = 1, . . . ,k. We shall estimate 
Pki {xi ,X2t . ■ . , Xki ; s), the other / — 1 factors are estimated in the same way. 

First assume that none of the variables xi, X2, . . . , Xki are within distance 
s from each other. Then one can clearly estimate pk^ {xi, X2, . . . , Xk^ ; s) from 
above as 
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PkAxi,X2, ■ ■ ■ ,Xk,;s) < / .../ p2kAxi^ ■ ■ ■ iXki,yi, ■ ■ -ykjdyi . . .dyk,- 

Jxi Jxk^ 

(48) 

Now, since p2ki{xi, . . . ,Xki,yi, ■ ■ - yki) = K[xi, . . . ,Xki,yi, ■ ■ ■ ,yfej, and 
K[xi, . . . , Xki,yi, ■ ■ ■ ,yki] is the determinant of a (2fci)-dimensional (non- 
negative definite) real symmetric matrix, we can estimate the determinant 
from above by the product of the determinants of the 2x2 diagonal blocks 

fci 

K[xi, . . .,xki,yi, ■ ■ ■,yki] = K[xi,yi, . . .,Xk^,yk^] < Y[K[xi,yi\. (49) 

i=l 

The bound H49|) follows from the Fischer inequality which we state below as 
Lemma 2. 

Lemma 2. Let M = ( ^] be a block matrix, let A and C be n x n and, 



B* C ^ 

respectively, m x m non-negative definite matrices, and B be a m x n matrix. 
Then 



det ^ ^ ^ - '^^^ ^ '^'^^ 

Proof 

To prove Lemma 2, it is enough to reduce it to the obvious case M = 
^B* by appropriate rotations and dilations in C" and C" (see e.g. 

CHI). 

As was shown above (see calculations around formula ^60\ 

rXi+s 

/ K[xi,yi]dyi < const{s)^, (51) 

J Xi 

which then implies that 

Pk, {xi,X2,...,Xk^) < const''' (S)3'=i . (52) 

If none of the variables are within s from each other in all factors in (|47|l . 
We infer from that 

Pki (xi , X2, . . . , Xfe, ; s)pfc2 (a^fci+i, ■ ■ • , Xk^+k2 ;s) x ... 

Pk, (a;fei+...fc,_i+i, . . . , Xfe; s) < Constisf" = 0{L-''), (53) 

and the integral of the l.h.s. of (|53|l over [0, L]'' f] D'^ goes to zero as L ^ cxd, 
since ^^([0, Lf f] D") = O(i'^-i). 

If some of the variables in pk^ {xi, . . . , Xk^ ) are within the distance 5 from 
one another, the analysis is quite similar. Let us assume, for example, that 
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Xi < X2 < ■ ■ ■ < a^fei , and that Xi < x^+i < Xi + s, i ~ 1, . . . ,p, and that 
the rest of the variables Xp+i, . . . , x^-^ are not within the distance s from each 
other. Then 

PkAxi,x2, ■ ■ ■ ,xki;s) < / ... / p2ki~p{xi, ■ ■ ■ ,Xki,yp+i, ■ ■ ■yki)dyp+i . . .dyk^. 

(54) 

One can write 

P2ki-pixi, ■ ■ .,Xki,yp+i,. -.yki) = K[xi, . . . .Xk^.y^+i, ■ ■ • ,yfci] < K[xi,X2, ■ ■ . ,Xp+i,yp+i] 

kl 

X W K[x,,y,]. (55) 

i=p+2 

As before, 

/ K[xi,yi]dyi < const{s) , i = p + 1, . . . , ki. (56) 

J Xi 

As for the term K[xi, . . . , a^p+i, y^+i], one can substract the first column from 
all other columns, and obtain 

K[xi, . . . ,Xp+i,yp+i] < Const kisf'^^, (57) 

since \g{x) — g{y)\ = 0{s'^) ior < x,y < s (we used the fact that g'{0) — 0). 
Combining 15t)|) and (|57|l . and integrating over the t/'s we obtain 

PkAxi,X2, ■■■,Xk,;S)< Const{Sf'''-P+^ = 0{L-^'+i-i). (58) 
Note, however, that 

Vol{{xi, ...,Xk,):x,< x,+i < x,+rs, t^l,...,p} f|[0, if' = 0(L'=i-f ), 

(59) 

and the product of the right hand sides of and itS^ goes to zero. 

If there are several factors in (|53|) for which there are variables within 
distance s from each other, the analysis is very similar, and we leave the 
details to the reader. Combining all the estimate together, one concludes the 
integral of the l.h.s. of (|53|) over [0, L]"^ f] goes to zero as L ^ oo. This 
finishes the proof of Lemma. 



The result of Lemma 1 and formula (|18|l imply that limi^oo J2n=i 



+00 C„(L) 



X 



= as^{e^ — 1), where {C„(L)}^J^ is the sequence of the cumulants of the 
counting random variable N{L), where N{L) is the number of the points of the 
s-modified random pont field in the interval [0, L] . It follows from the definition 
of the s-modified random pont field that N{L) ~ Ni{L)+N2{L), where Ni{L) 
counts the number of particles of the original random point field that have 
exactly one neighbor within distance s to the right, and N2{L) counts the 
number of particles of the original random point field that have more than 
one neighbor within distance s to the right. We claim that the probability 
that N2{L) 7^ is going to zero as L cx). Specifically, we establish 
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Lemma 3. 

lim EN2{L) = 0, (60) 

L — ^oo 

where E denotes the mathematical expectation. 

Since N2{L) is a non- negative, integer- valued random variable, H6()|l implies 
that Pr(iV2(L) 7^ 0) &s L ^ 00. 

The proof of Lemma 3 is elementary. We use the estimate 

pL I px+s i-x+S \ 

EN2{L)<j^ U J P3ix,yi,y2)dyidy2jdx. (61) 

As before one can show that p3{x,yi,y2) — K[x,yi,y2] — 0(s'*), and thus 
EN2iL) < {SfL = 0(1). 

Theorem 1 is proven. Theorem 2 immediately follows from Theorem 1. 
Indeed, the event {i]{L) > s} is exactly the event that there are no nearest 
spacings smaller than sL~3 between the particles in [0, L]. 

Acknowledgment Research was supported in part by the NSF Grant 
DMS-0405864. 



References 

1. A.Borodin and G.Olshanski, Distribution on partitions, point processes, and 
the hypergeometric kernel, Commun. Math. Phys., 211, 335-358, (2000). 

2. O.Costin, J. Lebowitz, Gaussian fluctuation in random matrices, Phys. Rev. 
Lett, 75, 69-72, (1995). 

3. D.J.Daley, D.Vere- Jones, An Introduction to the Theory of Point Pro- 
cesses, Springer- Verlag, New York, 1988. 

4. P.Diaconis, S.N.Evans, Immanants and finite point processes, J. Combin. The- 
ory Ser. A, 91, No. 1-2, 305-321, (2000). 

5. H.-O. Georgii, H.J.Yoo, Conditional intensity and Gibbsianness of determinan- 
tal point processes, J. Stat. Phys., 118, no. 1-2, 55-84, (2005). 

6. K.Johansson, Discrete polynuclear growth and determinantal processes, Com- 
mun. Math. Phys., 242, no. 1-2, 277-329, (2003). 

7. K.Johansson, Determinantal processes with number variance saturation, Com- 
mun. Math. Phys., 252, no. 1-3, 111-148, (2004). 

8. J. B. Hough, M.Krishnapur, Y. Peres and B. Virag, Determinantal processes and 
independence, arXiv:math. PR/0503110 preprint. 

9. R. Lyons, Determinantal probability measures, Publ. Math. Inst. Hautes Etudes 
Set., 98, 167-212, (2003). 

10. E.Lytvynov, Fermion and boson random point processes as particle distribu- 
tions of infinite free Fermi and Bose gases of flnite density, preprint arXivimath- 
ph/0112006. 

11. J.L. Lebowitz and J.K.Percus, Long-range correlations in a closed system with 
applications to nonuniform fluids, Phys. Rev., 122, 1675-1691, (1961). 



16 Alexander Soshnikov 



12. R. Lyons, J.Steif, Stationary determinantal processes: phase multiplicity, 
Bcrnoillicity, entropy, and domination, Duke Math J., 120, No. 3, 515-575, 

(2003). 

13. O. Macchi, The coincidence approach to stochastic point processes. Adv. Appl. 
Prob. 7, 83-122, 1975. 

14. O. Macchi, The Fermion process-a model of stochastic point process with re- 
pulsive points. Taransactions of the Seventh Prague Conference on Information 
Theory, Statistical Decision Functions, Random Processes and of the Eighth 
European Meeting of Statisticians (Tech. Univ. Prague, Prague, 1974), Vol. A, 
391-398. Reidel, Dordrecht, 1977. 

15. Y.Peres, B. Virag, Zeros of the i.i.d. Gaussian power series: a conformally in- 
variant determinantal process, Acta Math., to appear. 

16. D. Ruelle, Statistical Mechanics, W.A.Benjamin, Inc, New York, 1969. 

17. T.Shirai, Y.Takahashi, Random point fields associated with certain Fredholm 
determinants. I. Fermion, Poisson and boson point processes, J. Fund. Anal., 
205, No.2, 414-463, (2003). 

18. T.Shirai, Y.Takahashi, Random point fields associated with certain Fredholm 
determinants. 11. Fermion shifts and their ergodic and Gibbs properties, Ann. 
Probab., 31, No. 3, 1533-1564, (2003). 

19. T.Shirai, H.J.Yoo, Glauber dynamics for fermion point processes. Nagoya Math. 
J., 168, 139-166, (2002) . 

20. A. Soshnikov, Determinantal random point fields, Russian Math. Surveys, 55, 
923-975, (2000). 

21. A. Soshnikov, Level spacings distribution for large random matrices: Gaussian 
fluctuations, Annals of Math., 148, 573-617, (1998). 

22. A. Soshnikov, Gaussian limit for determinantal random point fields, Ann. 
Probab., 30, No. 1, 171-187, (2002). 

23. A. Soshnikov, Janossy densities of coupled random matrices, Commun. Math. 
Phys. , 251, 447-471, (2004). 

24. H.J.Yoo, Gibbsianness of fermion random point fields, preprint (2004). 



